Interplay between Zeeman interaction and spin-orbit coupling in a two-dimensional 

semiconductor system 
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We analyse the interplay between Dresselhaus, Bychkov-Rashba, and Zeeman interactions in a 
two-dimensional semiconductor quantum system under the action of a magnetic field. When a 
vertical magnetic field is considered, we predict that the interplay results in an effective cyclotron 
frequency that depends on a spin-dependent contribution. For in-plane magnetic fields, we found 
that the interplay induces an anisotropic effective gyromagnetic factor that depends on the orienta- 
tion of the applied field as well as on the orientation of the electron momentum. 
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The effects produced by different spin-dependent in- 
teractions in semiconductor structures are currently in 
the forefront of experimental and theoretical efforts in 
mesoscopic physics. The explosive activity is motivated 
by the desire for a deep understanding of quantum co- 
herence phenomena. The other driving force is the hope 
that the spintronics research would provide novel, low- 
dissipative microelectronic devices [l|, |2| . 

Most measurements of spin effects in semiconductor 
microstructures are performed applying a magnetic held. 
In this case, two intrinsic spin-dependent interactions 
naturally appear: the well-known Zeeman interaction, 
which couples the electron spin with the applied held, 
and the spin-orbit coupling. The latter has been exten- 
sively studied startingfrom two-dimensional electron gas 
to quantum dots |3|,|j,|li|6|,|jl|8j], and constitutes the basis 
of several proposals for spin-based applications [UEJEI- 

In a semiconductor, the spin-orbit coupling originates 
as a relativistic effect caused by electric fields present in 
the material (cf 01). Such a system can be described by 
the Hamiltonian 



l~t = T + TLd + TLr + Tiz 



(1) 



Here T = (p%, + p^j /2m* is the kinetic energy with the 
effective mass m* of electrons in the conduction band of 
the sample and p x , p y are components of the momentum. 
The electric held caused by the bulk inversion asymme- 
try of the crystalline structure contributes to the Hamil- 
tonian, Eq.JlJ, with the Dresselhaus term, Hd E3- This 
term has, in general, a cubic dependence on the momen- 
tum of the carriers. For a narrow [0,0, 1] quantum well, 
it reduces to the 2D, linear momentum dependent term 
Wn = P (p x o~x ~Vy a y) /fi" Here, the cr's are the Pauli 
matrices, and 6 is the intensity of this interaction (cf 
Ref ill . The electric field caused by the structure inver- 
sion asymmetry (SIA) of the heterostructure generates 
the Bychkov-Rashba term, Hr 0] . Since in the asym- 
metric quantum wells the SIA comes from the vertical 



direction, the Bychkov-Rashba interaction Hr has the 
form: Hr = a (p v a x — p x o~v) /h, where a is the corre- 
sponding strength [l2l ]. 

The correlation between different spin-dependent in- 
teractions is a key point in the understanding of spin 
phenomena in semiconductors. In this work, we shall 
study a two-dimensional (2D) semiconductor system in 
a magnetic field, described by the Hamiltonian pay- 
ing special attention to the interplay between the Zeeman 
and spin-orbit interactions. We will show that depending 
on the orientation of the magnetic held and the strength 
of the spin-orbit coupling this interplay produces a dra- 
matic effect on the electronic spin structure of a 2D het- 
erostructure. 

We hrst consider the case of -B||[0, 0, 1], i.e., of the mag- 
netic field parallel to the growth direction (z) of most 
heterostructures. In this case the momentum p is re- 
placed by a generalized momentum P = p+ (e/c)A, 
where the vector potential A of the vertical magnetic 
field B = (0, 0, B) is taken in the symmetric gauge 
A = B/2(—y, x, 0). In the latter case the Zeeman term is 
given by: Ti.z = g* ^bBg z /2 = e z a z /2, where g* rep- 
resents the bulk effective gyromagnetic factor, that is 
assumed to be constant, \ir = eh/2m e c is the Bohr's 
magneton and B the intensity of the applied field. 

When a single source of spin-orbit coupling is con- 
sidered, the Rashba or the Dresselhaus one, it can be 
completely solved 0, 0, ^j- The Hamiltonian com- 
mutes with the total angular momentum operator ( J z — 
L z + S z ) when the Bychkov-Rashba interaction is only 
considered. On contrary, when the Dresselhaus term 
presents alone, the symmetry operator is given by L z — S z 
jlq . If both terms are considered, no symmetry is avail- 
able and a solution is unknown. Below we propose an 
approximate solution for the full problem, which is one 
of the main results of the paper. 

In the absence of any spin-orbit coupling, one obtains 



2 



Landau levels spaced in energy by the cyclotron fre- 
quency (lo c — eB/m*c) and having a spin splitting given 
by the Zeeman interaction. To include the effects of the 
spin-orbit coupling, we transform the Hamiltonian Eq.JQ) 
to the reference frame where it is diagonal in a spin space 
up to a second order 



corresponding orbital effect of the magnetic field. In this 

way, the last term in Eq.JSJ) is transformed to a kinetic- 
2 



like one 



(2) 



The explicit expressions of the operators involved in the 
unitary transformation read, 



O = P x cj x 



Py®y : Q 



and 
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a/h 



Py^x ~1~ Px&y 



ftWr 



huj r 



(3) 



(4) 



We expand the transformed Hamiltonian, Eq.J2J), up 
to the second order in the parameters of the transforma- 
tion. Third order terms can be neglected if the magnitude 
of the operator involved in the transformation is much 
smaller than unity. This requirement is satisfied for the 
spin-dependent terms that are dominated by the orbital 
effect of the magnetic field. Therefore, the validity of 
the present approach is roughly restricted to the regime 
where the energy scale concerning the spin-dependent in- 
teractions (e 2 , a, (3) is much smaller than the energy scale 
given by the cyclotron frequency (htOc)- Nevertheless, 
note that no restriction is imposed on the relative weight 
of the spin-orbit coupling and the Zeeman interaction. 

Under the above hypothesis, the resulting transformed 
Hamiltonian is diagonal in spin space 



w = 



I 



2m* 

e z 



m 

1 - T7T 



a 



a 2 



in 



H 2 \ hw r — e z 

2 



o 



hw c + e z 



(Pi + Pi 



2 h 2 \l-e z /fuv c l + e z /huj c 
a/3 e z 



2o/3e^ 



( h ^y-et ) ( P l + P y) a - The intensity of 
this term has a sixth order dependence on the param- 
eters of the spin-dependent interactions (ot,0, e z ). It is 
much smaller than the first kinetic terms in Eq.©. As 
a consequence, this term has no relevant effect on the 
spectrum and can be safely neglected. 

Under this approximation, one can see that the main 
effect of the spin-orbit coupling and its interplay with the 
Zeeman interaction is contained in the first term of the 
Hamiltonian JSJ. In the transformed reference frame, 
this term is composed of the initial generalized kinetic 
term plus a spin-dependent one. The third term is just a 
constant that redefines the origin of the energy, whereas 
the second one corresponds to the Zeeman interaction 
associated with the vertical field. 

Comparing with the usual expression for the kinetic en- 
ergy, one observes that the spin-orbit interaction cause a 
redefinition of the effective mass of the system. This ef- 
fective mass (m*- 1) depends now on the spin orientation 
of the electrons as well as on the parameters characteris- 
ing the different spin-dependent interactions. In this way 
we face with a Landau problem of free electron motion 
(with this effective mass) in the magnetic field. 

The spectrum associated with the Hamiltonian, Eq. JSJ 
without the last term, can be obtained straightforwardly 
and it reads 



e n i s = (2n + |/| + I + l)hu CtS + —s 



m 



a" 



il 2 



where 



fku c 



m 



l-e z /fruj c l + e z /hw t 



i3 2 



l-e z /hu) c l + e z /hiu c 



(7) 



(8) 



h 2 e 2 — {hw c )' 



-(P x Py+PyP x )a z 



and the corresponding eigenstates have their spin 
aligned in the vertical direction (a z ). The last term in 
Eq.lJSJl represents the contribution of the interplay be- 
tween the whole set of spin-dependent interactions: Dres- 
selhaus, Bychkov-Rashba and Zeeman. It is diagonal in 
spin, although it has a crossed dependence on the in- 
plane components of the generalized momentum. To illu- 
minate the relevance of this term to the spectrum, we ap- 
ply a second transformation to the Hamiltonian, Eq. J5j), 



af3 



H 2 facie* - (fac) 2 ) 



(6) 



and keep all contributions, up to the second order. 
This approximation is well justified within our approach, 
since the magnitude of the term is much smaller than the 



and s = ±1 denotes the different spin states. In Eq.Q, 
(5) the expression '(2n + |/| + I + 1)' is usually replaced by 
'(2N +1)' corresponding to the principal Landau level 
index. The second term is the usual Zeeman contribu- 
tion and the last one represents the constant shift. Note, 
that the correction introduced by the spin-orbit coupling 
in the spectrum, Eq.©, depends quadratically on the in- 
tensities of the different spin-orbit terms, in contrast with 
the situation at zero magnetic field, where the effect in 
the spectrum depends linearly on these intensities. For a 
quantum dot, when only Zeeman and Rashba terms are 
considered |I9|, the interplay between the above terms 
modifies the spin-orbit intensity, while the kinetic term 
manifests a weak dependence on the spit-orbit coupling. 
For g* — (Zeeman term is absent) in a quantum dot, 
the kinetic term and the spin-orbit coupling remain with- 
out modifications. 

When the Zeeman term is absent in the 2D system, 
the resulting spectrum shows, however, an effective spin- 
dependent cyclotron frequency, U! CjS — uj c — 2m* (a 2 — 
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9* 


m*/m e 


£ z /htO c 


GaAs 


-0.44 


0.067 


-0.015 


InAs 


-14.9 


0.023 


-0.169 


InSb 


-51.6 


0.014 


-0.355 


AlAs 


1.52 


0.15 


0.112 


Ino.53Gao.47As 


-4.38 


0.038 


-0.082 



TABLE I: Characteristic values of the parameters involved 
in the ratio between the Zeeman energy and the cyclotron 
frequency for different semiconductor materials. 



f3 2 )s/H 3 . Materials with the effective g* — can be cre- 
ated by changing the concentration of one of the compo- 
nents of a sample (see, for example, poj V Comparing the 
latter expression with that given by Eq. @, it can be 
seen that the different spin-orbit terms (Bychkov-Rashba 
and Dresselhaus) interact in a different way with the Zee- 
man one. In both cases, the effect is quantified by the 
ratio between the Zeeman energy and the conventional 
cyclotron frequency (=e z /hui c ). Note, that this ratio is 
independent on the strength of the magnetic field, since 
it can be expressed in terms of the gyromagnetic factor 
and the effective mass 



off 



(9) 



Therefore, the interplay acts, in an effective way, renor- 
malising the intensities corresponding to the different 
spin-orbit terms. The character of this renormalisation 
depends as well as on the sign of the gyromagnetic factor 
of the semiconductor material. If the sign of the gy- 
romagnetic effective factor is negative then the effective 
strength of the Bychkov-Rashba term is reduced, whereas 
the corresponding to the Dresselhaus term is enhanced. 
On contrary, if the sign of '<?*' is positive, the effective 
Bychkov-Rashba term is enhanced and the Dresselhaus 
one tends to be suppressed. Thus, our analytical results 
provide a transparent explanation for numerical results 
obtained for 2D electron gas with Rashba spin-orbit cou- 
pling only pH and predict various new phenomena re- 
lated to the interplay between both spin-orbit and Zee- 
man terms. 

In table [I] there are listed the typical values of the ef- 
fective gyromagnetic factor and the effective mass corre- 
sponding to different semiconductor materials. A wide 
gap material such as GaAs, characterized by a relatively 
weak spin-orbit interaction, exhibits a small correction of 
the spin-orbit effect of about 1.5%. However, in the case 
of narrow gap semiconductors, which are characterized 
by strong spin-orbit intensities and g-factors, the inter- 
play between both interactions yields a sizeable correc- 
tion of the effect corresponding to the different spin-orbit 
terms of about 17% for InAs and 35% for InSb. 

If the magnetic field is parallel to the plane of an elec- 
tron motion, the situation changes substantially. In 2D 
systems, the orbital effect of the applied in-plane mag- 
netic field is frozen due to the strong confinement in 



the vertical direction. The only effect of the field ap- 
pears through the Zeeman interaction. In this case, the 
Hamiltonian, Eq.Q, preserves the translational invari- 
ance and the momentum is conserved. The Zeeman term 
now reads 



TLz = y (cos 8a x 



sin 8a y ) , 



(10) 



where the angle 8 represents the azimuthal orientation of 
the magnetic field in the x — y plane. The corresponding 
spectrum is given by the following analytical expression: 



A - 
B 



h 2 k 2 

2m* 
■ ak v - 

- ak x - 



f s[A 2 + B 2 ] 1/2 

(3k x + e z cos 8/2 
- [3k y — e z sin 8/2 



(11) 



where 'k x ' and 'k v ' are the components of the in-plane 
momentum. This is a contrast to a quantum dot, where 
the approximate analytical solutions are available only 
for Hz = and for Hn,T> <^H Z .18]. 

To elucidate the implications of the above expression, 
let us analyse some illustrative limit cases. First, we will 
consider the case when a single source of spin-orbit cou- 
pling is present. In the absence of the magnetic field, the 
spectrum is defined by — h 2 k 2 /2m* ± o^fe (a, = a, (3). 
The spin splitting 'ajfc' is isotropic, i.e., it doesn't depend 
on the orientation of the momentum. However, this sym- 
metry disappears when the Zeeman term is included. In 
the limit (3 = the spin splitting can be written as 



(^ak + y ^ + ae z k(sin((f> — 8) — 1) 



-1 1/2 



(12) 



where the angle '0' represents the azimuthal orientation 
of the momentum in the x-y plane. Note that, when 
the momentum is perpendicular to the magnetic field 
(4> = tt/2 + 8), the spin-orbit coupling and the Zeeman 
terms are decoupled in the spectrum. The same happens 
for the Dresselhaus interaction, but, in this case, when 
the momentum and the applied field satisfy the condition 
(9 = —(f>). In general, however, the interplay between 
the Zeeman interaction and the spin-orbit coupling is re- 
flected by means of Eq. I|12|) . 

Let us consider the limit of small Zeeman energy com- 
pared to 'ah', which is usually satisfied near the Fermi 
level of typical semiconductor heterostructures. Expand- 
ing Eq. JT^J, we obtain 



h 2 k 2 



^ k x , k-u , s 



2m* 



+ s ak + g* eff { 



>¥)■ 



where 



9eff( 



9 



(13) 



(14) 



is an anisotropic effective gyromagnetic factor originated 
due to the interplay between the two spin-dependent in- 
teractions. Thus, at B 7^ 0, the splitting partners, de- 
termined by the term sak in Eq. i|13|) . are enriched by 
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the second term that is linear dependent on the mag- 
netic field. Due to a random orientation of an electron 
momentum in plane, the borders of such sub-bands are 
restricted within the range [— <7*/xs_B/2, +g* HbB /2], cor- 
responding to the different momentum orientations , <j>\ 
instead of the discrete spin states ±g* /j,bB /2 associated 
with an isotropic g-factor. Depending on the relative ori- 
entation between momentum and the applied field, the 
value of the effective g-factor can change its sign with re- 
spect to the initial bulk value. In fact, in experiments 
with a tilted magnetic field [2^ it was observed that 
the effective g-factor depends on the tilted angle in the 
growth direction (z). In this case the increase/decrease 
of the g-factor may be influenced by the orbital motion 
alone, since the tilting of the magnetic field affects dif- 
ferent components of the orbital momentum. In contrast 
to this experiment, we provide a transparent mechanism 
of the anisotropy of the effective g-factor in plane due 
to the interplay between magnetic field and strong spin- 
orbit coupling. This phenomenon may serve as a guide- 
line to construct a novel class of read/write devices based 
on nanostructures with the strong spin-orbit interaction. 

In the opposite limit e z 3> cnjfe, the spectrum can be 
approximated by 

h 2 k 2 ( , s Q*UbB\ , s 
s kxMy . s ^—+ S [aks^-0) + y -^— J. (15) 

Consequently, in this limit, there is an angular de- 
pendence of the orbital contribution to the spin split- 
ting, which can be attributed to an anisotropic effective 
strength of the Bychkov-Rashba term. 

Another interesting limit is that the intensities of both 
spin-orbit terms are equal (|a| = Assuming a = —fi, 
the spin splitting reads 

r 2 I 1 / 2 

s (^2ak yx + + V2ae z k yx (cos (o - ^ - l) 

(16) 



where k yx — k y — k x . When the magnetic field is ori- 
ented at 9 = 7r/4, the spin-orbit and the Zeeman terms 
don't couple in the spectrum, independently on the mo- 
mentum's orientation. In this case the Hamiltonian has a 
symmetry in a spin space and, consequently, orbital mo- 
tion and spin are effectively decoupled and the spin-orbit 
interaction is effectively suppressed. Thus, we gene ralize 
the result found at B = by Schliemann at el |ll| for a 
nonzero magnetic field. 

In summary, we have analytically studied the effect 
of an applied magnetic field on the interplay between 
the Zeeman, Bychkov-Rashba and Dresselhaus interac- 
tions in a 2D electron system. When the magnetic field 
is perpendicular to the plane, we found that the spin- 
orbit interaction induces a renormalization of the effec- 
tive mass that becomes dependent on the spin orienta- 
tion of the electrons. This spin-dependent contribution is 
noticeably influenced by the interplay between the Zee- 
man interaction and the spin-orbit coupling, which mod- 
ifies in a different way the effective strength correspond- 
ing to each spin-orbit term. The spin-orbit coupling ef- 
fectively transforms to a spin-dependent cyclotron fre- 
quency (w CjS = eB z /m*c) that imposes different spatial 
scale of an electron motion. Indeed, the spin separation 
of cyclotron motion has been observed recently in j23| but 
in case of 2D hole gas, which exhibits a stronger spin-orbit 
interaction. The magnitude of the interplay between 
both spin-dependent interactions has a sizeable value for 
a narrow gap semiconductors with large g- factors. For 
the in-plane magnetic field, the interplay induces the 
anisotropic effective gyromagnetic factor that would lead 
to the observation of a continuum of Zeeman sub-levels 
in the spin structure of the spectrum. When the Zee- 
man interaction dominates over the spin-orbit coupling, 
an anisotropic effective spin-orbit intensity emerges. 
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